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NUCLEAR MAGNETIC RESONANCE METHOD OF DETECTING AND 
MONITORING THE FLOCCULATION KINETICS OF HEAVY FRACTIONS OF A 
COMPLEX FLUID 

BACKGROUND OF THE INVENTION 

Field of the Invention 

The invention relates to a nuclear magnetic resonance method of detecting and 
monitoring the flocculation kinetics of high molecular weight fractions of a complex 
fluid. 

Description of the Prior Art 

Flocculation and deposition processes pose considerable problems in the petroleum 
industry. In particular for heavy oils, components of very high molar mass (asphaltenes, 
resins) are often the cause of such processes which may appear in porous media during 
production as well as during transportation. Flocculation is the formation of molecular 
aggregates of micronic size leading to sedimentation or deposition that can considerably 
modify the fluid flow, either by reduction of the section of flow or by viscosity increase. 
Furthermore, the intrinsic charge of some components (for example asphaltenes) 
generates a high tendency to cling to the charged surfaces. 

The thermodynamic parameters which govern the flocculation processes are 
numerous (composition, pressure, temperature) and the complexity of the molecular 
structures involved make prediction and modelling very uncertain. Similarly, certain 
recovery methods (C02 injection, acidizing) may modify the fluid equilibria and bring 
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about these processes. It is thus necessary to carry out measurements but the available 
techniques do not allow the first stages of the flocculation process to be observed and 
pose considerable implementation problems as regards pressure and temperature or in- 
situ problems in oil wells. 

5 Many petroleum crudes, notably those referred to as asphaltene crudes, are liquid 

hydrocarbon fluids which contain more or less large amounts of heavy fractions in the 
dissolved state and/or in the stable colloidal state in the pressure and temperature 
conditions to which the fluids are subjected. When these pressure and/or temperature 
conditions vary, notably when the pressure decreases, the heavy fractions contained in 

10 these fluids tend to flocculate and to settle in the formation in the neighbourhood of 
wells, in wells and in production and transfer facilities intended for the fluids. Thus, 
when a hydrocarbon reservoir containing heavy fractions is developed, generally before 
the bubble point is reached, the stability of these fractions decreases. When the 
saturation threshold is reached, the heavy fractions flocculate and settle, which can 

15 cause clogging of the porous media and formation of plugs likely to severely damage 
production wells and surface installations. 

For oil producers whose task is to extract and convey, through production wells and 
pipe networks, liquid hydrocarbon fluids of petroleum crudes containing heavy 
fractions, for example asphaltene crudes, from production fields, it is therefore 
20 important to have precise knowledge of the pressure thresholds below which the heavy 
fractions settle, so as to carry out production and transfer of the fluids under pressure 
and temperature conditions preventing deposition of the heavy fractions in installations 
or to provide a suitable treatment. 
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Various methods of determining the deposition threshold of heavy fractions, 
notably asphaltenes, contained in liquid hydrocarbon fluids of petroleum crudes are 
known. These methods are most often optical light transmission or diffusion methods, 
conductimetric methods or viscosimetric methods. 

5 A method known as the "spot test" deposits a small quantity of a mixture on a filter 

paper and observes the spot that forms. The flocculation aggregates that form in a 
mixture diffuse less readily than the surrounding liquid. Thus, if the spot is not uniform, 
it is an indication that it contains flocculating particles. 

The aforementioned methods use detection of the variation of a physical quantity, 
10 for example an absorption coefficient or absorbance of light rays in the visible range or 
in the infrared range, electrical conductivity or viscosity, which results from the change 
in the structure of the fluid following flocculation and deposition of heavy fractions. 

A major drawback of such methods is that they are not very selective insofar as it is 
not always easy to relate the variation of the physical quantity measured to the 
15 flocculation and deposition of heavy fractions, and these methods are not always 
sensitive to the deposition of a small amount of such fractions. Some methods, such as 
absorbance measurement in the infrared range, are very sensitive but difficult to 
implement under reservoir conditions. 

Furthermore, since these methods are often used in the laboratory, the question 
20 which has to be considered is the representativity of the samples on which the physical 
quantity measurements are carried out. In fact, for a sample to be representative of the 
sampled fluid, it is necessary to maintain this sample under the pressure and 
temperature conditions that prevail for the sampled fluid, for example reservoir fluid, 
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throughout the sampling, sample transport and storage operations preceding the 
measurements, 

French patent 2,818,753 filed by the assignee describes a method of determining the 
deposition threshold of heavy fractions contained, in the dissolved state and/or in the 

5 stable colloidal state, in a liquid hydrocarbon fluid. The invention provides a method of 
determining the deposition threshold of heavy fractions, notably asphaltenes, contained 
in the dissolved state and/or in the stable colloidal state in a liquid hydrocarbon fluid, 
using the formation of an increasingly high pressure drop linked with the flow, at an 
increasing flow rate, of a sample of the fluid through a capillary passage. The fluid 

10 sample being in the dissolved state and/or in the stable colloidal state, at the inlet of a 
capillary passage a pressure drop which is at least equal to the difference between the 
pressure of the fluid sample and the bubble-point pressure of said sample is generated 
between the inlet and the outlet. A significant shift in the variation as a function of time 
of AP (difference between the pressure of the fluid at the capillary inlet and the pressure 

15 at the outlet) and of a quantity D representative of the flow of liquid flowing through the 
capillary passage is detected, which allows to characterize the deposition threshold of 
the heavy fractions of the fluid. 

Besides, it is well-known that NMR devices can be used notably to measure certain 
physical characteristics of fluid mixtures such as hydrocarbons, notably the viscosity or 
20 the gas/oil ratio (GOR). The viscosity of the mixture and its GOR coefficient are 
obtained from the NMR measurement of a diffusion coefficient D and from the 
measurement of the longitudinal Tl or transverse T2 relaxation time. Such an 
application is for example described in WO-0, 142,8 17 or United States patent 
5,696,448, or in the following document : 
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- Prammer, M.G. et al., (2001);"the Downhole NMR Fluid Analyser"; SPWLA 42 nd 
Annual Logging Symposium. 

It can be noted that these measurements are generally performed in a time interval 
after excitation which is unsuited to detection and interpretation of flocculation 
5 phenomena. 

It is also well-known to use NMR devices to detect and monitor a very different 
phenomenon which is the crystallization of particles in fluids. 

SUMMARY OF THE INVENTION 

The inventors have observed that NMR type methods applied to the detection of 
10 solid particles can also surprisingly apply to non-solid particles of high molecular 
weight under slow rotation which progressively aggregate, and are checked so that their 
flocculation rate can be determined by means of this type of method. 

The application of the method according to the invention is notably to monitor the 
flocculation kinetics of generally asphaltene polar fractions which are contained in the 
15 dissolved state and/or in the stable colloidal state in a liquid hydrocarbon fluid. 

The nuclear magnetic resonance method of detecting and monitoring the 
flocculation kinetics of non-solid high molecular weight aggregates of a complex fluid 
comprises applying to the fluid a first static polarization magnetic field, then at least a 
second oscillating pulsed magnetic field perpendicular to the first one, created by coils 
20 connected to an excitation generator for nuclear magnetic resonance of the nuclei being 
considered and acquisition of the relaxation signals of the nuclei in the fluid. 



The method comprises detecting, regarding the relaxation signals, a first part 
representative of the relaxation of these aggregates in the fluid and a second part 
representative of the relaxation of the liquid fraction of the fluid, and determining the 
flocculation rate (Tf) of the aggregates by comparison of the values extrapolated at the 
start of the acquisition times of the first part and of the second part respectively. 

The flocculation rate can be determined by means of the relation : 

Tf = (M x (t = 0) - M xl (t = 0))/M x (t = 0) 

where M x (t = 0) and M xl (t = 0) are the values extrapolated at the start of the acquisition 
times of the first part and of the second part respectively. 

The flocculation threshold of the fluid can for example be obtained by modelling 
the relaxation signals actually obtained by means of a combination of exponential 
functions depending on an adjustment parameter and the threshold corresponding to a 
maximum value of said adjustment parameter is sought. 

According to an implementation mode, the method comprises applying to the fluid 
a sequence of two 90° pulses referred to as pseudo-solid echoes in which a 180° 
magnetization focussing pulse is inserted, between two successive applications of the 
90° pulses, with time intervals t/2 between the different pulses, and measuring the 
maximum amplitude of the relaxation signals in the neighbourhood of time t=2x for 
different values of t in the sequence. 

The method affords many advantages. It allows continuous monitoring, useful for 
analysis of the flocculation as a function of the chemical composition and of the solvent. 
All of the volume isanalyized. Unlike optical methods, even the non-transparent 



samples can beanalyized. Characterization of the kinetics is easy. It is also easy to 
evaluate the aggregate ratio, which can lead to an approximate estimation of the 
molecular weight. 

BRIEF DESCRIPTION OF THE DRAWINGS 

5 Other features and advantages of the method according to the invention will be 

clear from reading the description hereafter of non-limitative embodiment examples, 
with reference to the accompanying drawings wherein : 

- Fig. 1 shows the transverse magnetization decrease measured on a flocculated mixture, 

- Fig. 2 shows the decrease curves of the transverse magnetization measured for an 
10 initially non-flocculated, then flocculated mixture, 

- Fig. 3 diagrammatically shows a NMR type analysis device, 

- Fig. 4 shows the relaxation signal of an asphaltene/toluene/heptane solution, and 

- Fig. 5 shows the evolution of the coefficient A of a function modelling the relaxation 
signal of Fig. 4, which allows to locate the value of the flocculation threshold. 

1 5 DETAILED DESCRIPTION OF THE INVENTION 

NMR analysis technique essentially applies to an object to be tested a first static 
polarization magnetic field Bo intended to align the nuclei of the initially randomly 
oriented hydrogen protons, in the direction of the field, then a second oscillating pulsed 
magnetic field at the Larmor frequency, perpendicular to the first one, created by coils 
20 excited by a control signal to carry out a nuclear magnetic resonance experiment. When 
this pulsed field stops, the return of the nuclei to their initial state or relaxation 
generates electromagnetic signals (echoes) that are detected and analysed. The presence 
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of such and such substance and some of its physical parameters are determined from the 
amplitude characteristics of these signals. 

In cases where the analysis concerns a fluid mixture, line 1 in which it circulates (or 
possibly the vessel that contains it) is conventionally passed (Fig.3) first in the air gap 
5 of a permanent magnet 2 in which a static polarized field B 0 (z) prevails in a direction z, 
then in a coil 3 connected on the one hand to a generator 4 delivering on demand a 
current pulse at a frequency in the radio frequency range creating an oscillating pulsed 
transverse field B,(t) in a direction x perpendicular to direction z, and on the other hand 
to a relaxation signal detection circuit 5. 

10 Within the scope of the method according to the invention, analysis of the 

relaxation of the transverse magnetization M x (t) of the protons of a liquid considered, 
observed at short times, is applied to detection of the flocculation, to determines the 
proportion of coarse structures formed in suspension in a liquid and to characterize the 
kinetics governing this formation process. 

15 The procedure can be carried out conventionally using a sequence referred to as 

Hahn sequence or a sequence referred to as pseudo-solid echo sequence. After 
polarization by the static magnetic field Bo(z) of the order of 0.5 T for example, a pulse 
field B^t) having for example a tilt front and a receiver desaturation time below some 
microseconds is applied by means of coils 3. A sequence of at least two 90° pulses 

20 separated by a time interval t is applied, with an inserted 180° pulse focusing the 
magnetization and thus being free from the diamagnetic heterogeneities of the material. 
The pulse application times are separated by time intervals x/2. The sequence can be 
schematized in time as follows : 
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(90°)x - x/2 - (180°)y - x/2 - (90°)x - echo acquisition. 

By analyizing the echoes by means of a device 6, their maximum amplitude M x (t) in 
the vicinity of time t=2x is determined for different values of t in the above sequence. 

For a liquid in which certain fractions have flocculated, the relaxation signal has 
two clearly distinct parts corresponding to the appearance of the coarse structures in the 
liquid : a first part decreasing very rapidly within about twenty microseconds and whose 
size increases regularly with the formation of the flocculated structures, and a second 
part decreasing very slowly within several hundred milliseconds and corresponding to 
the liquid part of the fluid (Fig.2). The flocculated structures produce fast decrease of 
the magnetization because they can be considered to be congealed (nearly solid). They 
are thus clearly distinguished from the magnetic properties of the surrounding liquid. 
Besides, the flocculated structures are also subjected to the Brownian movement 
generated by the thermal agitation of the liquid, but this random movement is too long 
to induce a magnetic relaxation process. This approach applies when the average size of 
the structures exceeds some nanometers. 

The flocculated aggregate rate Tf can be obtained by extrapolating the signal of the 
liquid part M xl at the time t=0. 

A graphic method can be used when, as in the example of Fig. 1 , the distribution of 
the points relative to the liquid part of the fluid analysed is substantially linear. It is also 
20 possible to seek for example a polynomial of sufficient degree or, if necessary, a sum of 
exponential functions modelling the distribution of the points of the liquid part, and to 
deduce therefrom the extrapolated value at the time t=0. 
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Rate Tf is obtained from this value M xl (t=0) and from the corresponding value of 
the first part of the distribution close to the start of the times, similarly extrapolated at 
the time t=0, i.e. M x (t=0). For example, Tf can be defined by the relation : 

Tf = [M x (t=0) - M xl (t=0)] / M x (t=0) (1) 

5 where M xl = M x (t>tc), tc being the time separating the two decrease regions of M x (t). It 
is implicitly assumed that the total mass analysed does not change. 

Detection of the first part of the relaxation distribution (Fig.l) implies that the 
relaxation signals are acquired practically from the start of the times, that is within an 
interval greatly below 1 ms, which is not the case with the conventional methods using 
10 NMR analysis devices. 

Observation of the flocculation kinetics as a function of time is illustrated in Fig. 2 
for a viscous liquid (100 centipoise). At t = 5 days, there is no flocculation because the 
signal shows no fast decrease. At t = 29 days, certain flocculated structures appear and 
become more numerous at t = 36 days. The acquisition time for a curve M x (t) being 
15 some minutes, faster kinetics can readily be observed. 

The method application example below, which relates to asphaltene/toluene/heptane 
solutions, clearly shows a flocculation threshold detection mode. 

It is well-known that, when the proportion of heptane increases, this type of solution 
flocculates. Fig. 4 shows an example of transverse magnetization signals (Mx(t)) 
20 obtained with these solutions. Parts PI and P2 described above are connected by an 
intermediate part P3. In general, the signals are correctly described by a function which 
is a sum of exponentials. In the following function : 
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Mx(t) = exp[-(t/T)°] + A exp(-t/T 2S ) (2) 

the second term A exp(-t/T 2S ) is characteristic of the relaxation time of the 
heptane/toluene mixture. Equation 2 is useful for separating the high molecular weight 
part characterized by very short relaxation times (order of magnitude x) from that of the 
5 heptane/toluene mixture (order of magnitude T 2S ). The first term describes part PI and 
intermediate part P3. For the mixture considered, the high molecular weight structures 
are visible before flocculation. 

The value of A is obtained by determining the coefficients of relation 2 so that the 
modelled curve coincides with the curve of Fig. 4. 

10 The aggregation rate of this mixture is calculated by means of the following 

relation : 

Tf=(l-A)/A (3) 

This way of calculating Tf is equivalent to that obtained when applying relation 1 . 
Monitoring of the evolution of coefficient A or, in an equivalent manner, of Tf, shows 
15 that the flocculation threshold corresponds to the maximum of this coefficient, in 
agreement with standard techniques. 

Fig. 5 shows the evolution of coefficient A as a function of the proportion T/H of 
toluene/heptane when the asphaltene concentration is 10.7, 8.8 and 7.6 % respectively. 
It can be seen that the flocculation threshold at the heptane-toluene mass proportion 35- 
20 65 % shown by the arrow corresponds to the maximum of coefficient A. 
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The NMR measuring device used for implementing the method can be miniaturized 
in order to be installed in an oil well in combination with a formation-test tool of a well- 
known type in any other surface installation to monitor the flocculation phenomena. 
Whatever the analysis conditions, about 50 mg. substance are sufficient in practice. 



